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Abstract 

We attempt to understand the fate of spacelike gravitational singular- 
ities in string theory via the quantum stress tensor for string matter in 
a fixed background. We first approximate the singularity with a homo- 
geneous anisotropic background and review the minisuperspace equations 
describing the evolution of the scale factors and the dilaton. We then re- 
view and discuss the behavior of large strings in such models. In a simple 
model which expands isotropically for a finite period of time we compute 
the number density of strings produced by quantum pair production and 
find that this number, and thus the stress tensor, becomes infinite when 
the Hubble volume of the expansion exceeds the string scale, in a manner 
reminiscent of the Hagedorn transition. Based on this calculation we argue 
that either the region near the singularity undergoes a phase transition 
when the density reaches the order of a string mass per string volume, or 
that the backreaction of the produced string matter dramatically modifies 
the geometry. 

1 Introduction 

Infinities and singularities in physical theories are signals that the theory in 
question is inconsistent or incomplete, and force physicists to uncover deeper 
structures. Much of twentieth century physics has arisen from such singulari- 
ties; two notable examples are the ultraviolet catastrophes of classical thermo- 
dynamics which helped drive the development of quantum mechanics, and the 
infinities in field theory loop calculations which helped drive the development 
and understanding of the rcnormalization group. 

Among the theories which we currently use to describe the world, general 
relativity remains doubly afflicted: at the quantum level it is nonrenormalizable, 
and at the classical level it contains singularities, as demonstrated by the well- 
known singularity theorems of Hawking and Penrose (Penrose 1965, Hawking 
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1967, Hawking and Penrose 1970). With both of these singularities we are faced 
with the problem of predictability. Nonrenormalizability tells us that there must 
be a very different description of the world at the Planck scale, but gives no clue 
as to what that description might be. Classical spacetime singularities defined 
via geodesic incompleteness tell us that in finite proper time, we will reach 
regions of spacetime beyond which we cannot evolve the dynamical equations of 
the theory. Of course, if we attempt to cure this classical problem by quantizing 
the theory, we must contend with the nonrenormalizability of the theory. 

One subject of active research which may force physicists to face these two 
problems directly is the apparent evaporation of black holes via quantum particle 
production (Hawking 1974). It is our belief that to understand the outcome 
of the Hawking process we must eventually understand the physics near the 
singularity. In particular this seems to be the primary lesson taught by the study 
of two-dimensional dilatonic black holes; at some point the singularity and the 
apparent horizon coincide and one needs some prescription for dealing with the 
singularity at this point in order to understand the final state of the black hole 
(see Strominger 1995 and references therein). To date there have been several 
proposals for the what endpoint the of black hole evaporation process might be|^ 
(for critical reviews and references see Banks 1994 and Strominger 1995); many 
of these proposals rely on low-energy physics up to the point that the black hole 
reaches the Planck mass, and then propose plausible consistent completions of 
the evaporation process, without recourse to a Planck-scale theory. Nonetheless, 
it remains to be shown explicitly that any of these proposals are the result of a 
complete theory of quantum gravity. 

Thus, if we claim to understand quantum gravity, we must understand and 
remove the singularities that seem to plague it. This is hardly a new enterprise, 
and while we are not capable of producing a full history of previous attempts 
to do so, one may divide the attempts into two categories: 

• Evasion of some tenet of the singularity theorems 

• Searching for a breakdown of the notion of classical spacetime in a region 
of strong curvature, via quantum mechanics or string theory 

The most natural way to avoid the singularity theorems is to violate the condi- 
tion 

7^a^,^C' > 0, (1) 

(where TZab is the Ricci tensor and ^ is a timclikc or null vector field depending on 
which theorem one appeals to); the other conditions involve global hyperbolicity 
of the spacetime (although this condition is non-essential), the lack of closed 
timelike loops, and the existence of trapped surfaces (which occurs given a 

^ "The number of very good physicists who have expressed fairly definitive opinions about 
the resolution of the Hawking puzzle is smaller than the number of definitive opinions they 
have expressed." (Banks 1994) 
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sufficient concentration of matter (Schoen and Yau 1983)), and are physically 
too attractive to casually dismiss. The singularity theorems are based on feeding 
Equation (|l|) back into the Raychaudhuri equation; in general relativity one can 
then use Einstein's equations to relate Equation to the strong or the null 
energy conditions. To avoid the singularity theorems one must then either 
violate the appropriate energy conditions, or suitably modify Einstein gravity 
so that Equation (|^) is violated. In particular, quantum fields violate the energy 
conditions (Epstein et.al. 1965). 

In the long run simply adding terms to Einstein's equations or relaxing 
energy conditions may not seem entirely satisfactory: we know that as we ap- 
proach large field gradients close to the singularity, quantum geometry and 
(sub)Planckian physics come into play. Still, at this point the burden is on 
the theorist to explain what exactly quantum geometry is, and then to show 
that quantum mechanics combined with one's favorite choice of a Planck scale 
theory is not also singular, especially given the nonrenormalizability of Einstein 
gravity. 

In this paper we would like to examine gravitational singularities using quan- 
tum string theory. To do this we should redefine singularities in terms of ex- 
tended objects and quantum mechanics; we can do this with fairly natural ex- 
tensions of geodesic completeness. Classically, we require that the solutions to 
the string equations of motion never cause the string to leave the spacetime in 
finite proper worldsheet distance. Quantum mechanically, we require that the 
first-quantized wave operator be well-defined and Hermitian everywhere. With 
these definitions we can discuss three types of singularities: singular sources of 
classical fields (i.e. charged point particles), timelike singularities, and spacelike 
singularities occurring (in a suitably defined coordinate system) at some point 
in time. 

Singular classical sources are easily dealt with: they arise generally from 
considering particles as point-like objects. In a first-quantized theory quantum 
mechanics smears out the wave function of a test particle enough to avoid any 
pathology at the source, so long as the singularity in the potential is weaker 
than 1/r^ (which is the case for the known forces of nature.) In quantum field 
theory calculations, we get divergences from large loop momenta and from small 
impact parameters. These are results of not properly regulating the theory at 
short distances. 

Timelike singularities, of which the aforementioned singularities are a sub- 
set, are also not generically troublesome. It is strongly believed that string 
theory has a minimum distance and thus soft high-energy behavior, due to the 
extended nature of the fundamental objects. This has been shown to resolve 
many potential timelike singularities in the theory, most notably in amplitudes 
for high-energy scattering (Gross and Mende 1987, 1988), orbifold singulari- 
ties (Dixon et.al. 1985, 1986), and singularities arising from the degeneration 
of Kahler structures in Calabi-Yau manifolds (Aspinwall et.al. 1993 and 1994, 
Witten 1993). This good behavior arises from the apparent lack of any opera- 
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tive definition of a point in string theory; no point on a given string can carry a 
finite fraction of the string's energy, so that any momentum transfer in a stringy 
scattering process is shared by the entire string. We should note that we do not 
want to smooth out all timelike singularities. Horowitz and Myers (1995) have 
argued that any stable theory with Einstein gravity as a low-energy limit must 
have solutions with timclikc singularities arising from singular initial conditions; 
in particular, a solution which looks like a negative-mass Schwarzschild solution 
at large distances should be singular even in a full theory of quantum gravity 
(remember that the total mass is measured at infinity, where only the Einstein 
term in the effective action will contribute). If we cannot throw away such so- 
lutions the theory becomes unstable since there is no stable ground state, i.e. 
empty space decays by pair-producing these negative-energy lumps (Horowitz 
and Myers, 1995). 

Spacelike singularities pose a genuinely serious and confusing dilemma, as 

they occur entirely because of the universally attractive nature of gravity past 
the event horizon, where the causal structure forces classical matter to collapse 
completely. These facts seem to be independent of the ultraviolet nature of the 
theory. So, while the singularity is an ultraviolet object and perhaps will be 
resolved by physics similar to that which resolves timelike singularities, their 
inevitable nature and the fact that they arise from a rather soft attractive force 
begs for a careful examination. At any rate it is safe to say that we do not 
have as much control over these singularities in string theory as we have over 
the known timelike singularities, and so those of us who hope that string theory 
describes nature must still confront spacelike singularities. 

The first words uttered by string theorists when faced with singular back- 
grounds is that they are simply signs that the expansion of the spacetime ef- 
fective action in powers of icun./^stn where icmv. is the characteristic length 
over which the curvature changes, breaks down. Yet there are known exact 
conformal theories that possess spacelike singularities, the best known exam- 
ple being the 2-dimensional SL{2,n)/U{l) black hole, (Witten 1991). One 
might then appeal to string loop corrections to the /3-functions (Lovelace 1986, 
Fischler and Susskind 1986, Callan et.al. 1987). While we will argue below 
that string loop effects are part of the story, we should be careful. For example, 
there are gravitational plane- wave backgrounds with null singularities, for which 
string propagation is singular (Horowitz and Steif 1990a and 1990b, dc Vega and 
Sanchez 1992). These backgrounds correspond to exact conformal theories at 
tree level (Amati and Klimcik 1989, Horowitz and Steif 1990a) because all in- 
teresting counterterms that one might produce from the Riemann tensor, the 
dilaton, or the antisymmetric tensor field, will vanish; for the same reason loop 
corrections to tlic background will also vanish. However, it should be noted 
that these null singularities by definition do not arise from nonsingular initial 
conditions. We should also recall the arguments of Horowitz and Myers (1985) 
mentioned above: if corrections to the low energy effective action generically 
tame singular backgrounds, we should worry that the theory is unstable. 
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The aforementioned string theorist would next point out that corrections to 
the effective action are not the end of the story; for example, strings propagating 
in orbifold geometries are quantum mechanically well behaved at the orbifold 
points. The lesson here is that we should relax our definition of singularities 
somewhat, and define them as regions of the background for which the wave 
operator for string theory becomes undefined, thus properly including orbifold 
geometries and some degenerations of Kahler geometries as nonsingular back- 
grounds. This definition is a natural extension of geodesic completeness as it 
defines a singularity as the point at which predictability, here defined as the 
ability to evolve quantum wave functions, breaks down. This relaxed definition 
does not save us from spacelike singularities. For example, the Lq operator in 
the SL{2, TV)/U{1) black hole backgroimd becomes ill-defined at the singularity: 
string wave functions {i.e. tachyon vertex operators) blow up at the singularity 
(Dijkgraaf et.al. 1992). In order to support our claim that a' corrections and 
the properties of quantum wave functions in string theory do not in general save 
us from spacelike singularities, we will in the appendix review arguments which 
indicate that the structure of higher-dimensional black hole singularities is de- 
scribed by a perturbation of the SL{2,TZ)/U{1) black hole. Note also that the 
caveats of Horowitz and Myers (1995) apply to quantum mechanical smoothing 
as well. 

None of this is to say that there is no solution to the problem of spacelike 
singularities in string theory, or that the answer does not involve the short- 
distance properties of string theory; rather, we simply wish to emphasize that 
the problem has not been solved and that it requires some care. In the light of 
these arguments, this paper is an attempt to reach a primitive understanding 
of the quantum dynamics of string field theory in singular backgrounds, as a 
precursor to understanding the full quantum evolution of the background. We 
will make this attempt by first examining the spacetimc background near the 
singularity and noting its generic and essential features, and then by examining 
string propagation in singular backgrounds, and finally by trying to crudely 
understand string field theory in a simplified spacetime background for which 
the scale factor changes over a brief period of time. We will find that the density 
of produced strings in a given state is exponentially suppressed by the mass 
of that state times the inverse rate of expansion, p^^, so that when this rate 
becomes as large as (J^]., the density of states overwhelms the suppression and an 
infinite number of highly energetic long strings are produced. One may object 
that it is not legitimate to quantize string field theory in a classical background 
spacetime without explicitly taking backreaction into account from the start; 
however, since string theory has an as adjustable dimensionless parameter the 
ratio istr/^pi far from the black hole, we may make this parameter small and 
see if string theory resolves the singularity in the weakly coupled region where 
we have some hope of knowing what to do. Our result indicates a sort of 
paradox, of which we see two possible resolutions: either stringy backreaction 
must smooth the singularity, or there is a Hagedorn-like phase transition as the 
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density of produced strings becomes comparable to a string mass per unit string 
volume. In the former case we will suggest that the combination of an equation 
of state for string matter derived by Gaspcrini et.al. (1991a, b) and the effects 
of the dilaton (Tseytlin and Vafa 1992) potentially cause spacetime to become 
non-singular. Without thinking too much, one could imagine a few scenarios. 
For example, the spacetime could stabilize as expansion or contraction cease, 
with the dilaton kinetic energy balancing the produced energy density as the 
spacetime runs to strong coupling. Another possibility is that the singular region 
relaxes to some de-Sitter like phase, as suggested by Frolov et.al. (1990) (sec also 
Poisson and Israel 1988) for the case of particle field theory, based on a limiting- 
curvature hypothesis for quantum gravity. One of us (Martinec 1995) proposed 
a realization of this latter possibility by appealing to T-duality: this paper 
arose from that work, although we no longer feel that the duality properties 
of string theory are necessarily relevant for our purposes. Our arguments will 
be crude and will ignore many issues such as the conformal invariance of the 
background, but we will argue that we have captured the physics of the effects 
we are interested in and that what we have ignored will not invalidate our 
results. 

Since our proposals depend on the time-dependent nature of the singulari- 
ties we are interested in, they may not apply to timelike singularities; the sorts 
of timelike singularities which Horowitz and Myers (1995) discuss may remain 
singular as required for the stability of flat space. As a limiting case, the singu- 
lar null-fronted plane wave spacetimes described by Horowitz and Steif (1990) 
and de Vega and Sanchez (1992) will have no particle creation (and no string 
creation) due to the presence of a null Killing vector, so that in this borderline 
case the singularity is untamed by our scenarios. 

The remainder of the paper shall be organized as follows: in Section 2 we 
will review features of singular backgrounds in general relativity and string 
theory. In the case of string theory, we will examine the cosmological equations 
of Tseytlin and Vafa (1992) with respect to various phenomcnological energy- 
momentum tensors that one may imagine, and discuss the singularities that 
may arise in these equations. In Section 3 we will describe string propagation 
in expanding and/or contracting backgrounds, reviewing the work of Gasperini, 
Sanchez, and Veneziano (Sanchez and Veneziano 1990, Gasperini et.al. 1991a, b). 
We will use a similar approximation to examine the propagation of large strings 
in a temporarily expanding model. We will also discuss the equation of state 
for string matter in these models. In Section 4 we will discuss free quantum 
string field theory in curved backgrounds, and describe techniques for computing 
approximations to the rate of string production in such backgrounds. In Section 

5 we will apply these techniques to a specific isotropically expanding model and 
find the aforementioned singularity in the rate of string production. In section 

6 we will discuss how one might extend this calculation to an anisotropic model 
and what features one might look for; we will also discuss how the produced 
matter, given the equation of state discussed in section 3, might afii'ect the 
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singularity. Section 7 will contain our conclusions and further speculations. 
Appendix A will discuss how one might relate the SL{2,TZ) /U{1) black hole to 
the Schwarzschild solution. 



2 Classical backgrounds in general relativity and 
string theory 

Belinskii, Khalatnikov, and Lifshitz (Lifshitz and Khalatnikov 1963, Belinskii 
et.al. 1970, 1982) have argued that the spacetime metric close to a spacclike 
singularity can be written as 



UaTlp) dx'^dx^ , 



(2) 



where pi, I, m, and n are functions of the spatial coordinates only, and the pi 
satisfy the equations 

3 3 
i=l i=l 

which are identical to the restrictions on the spatially independent exponents 
Pi in a vacuum Kasner spacetime: 



ds^ = -dt^ + J2t^P' {dx' 



(4) 



Equation is found by expanding Einstein's equations in a power series in t, 
assuming a power law singularity in the metric, with the leading exponents as 
shown: one finds that contributions from the intrinsic curvature of the spatial 
slices are subleading as compared to their extrinsic curvature (though there 
is the requirement that, if pi is the negative exponent, ? • V x I = 0, which 
means that we may use the direction specified by Z as a global coordinate) . This 
means that the spacetime is "velocity-dominated" in this region. If we look on 
a reasonably small spatial slice, the metric should look like the Kasner metric. 
The form (^) is certainly what we find near the Schwarzschild singularity: the 
metric inside the horizon (r < r^) for a D-dimcnsional Schwarzschild black hole 
is 



ds^ = - 




-1 













dt^ + r^dn^ 



D-2 



D-3 

- ( - 1 dr' 



D~3 



dt^ 



^dn 



D-2 



(5) 
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If we write T = —[2,rh/{D — l)](r/r/i)^-^ ^)/'^^ and p — t, we can rewrite the 
most singular part of the metric as 

-2(D-3)/{D-l) 



ds^ = -dT^ 



[D- 1)T 



2rh 



[D- 1)T 



2rh 



A/{D-l) 

JO 

D- 



(6) 

Note that in four dimensions the relations (^) are satisfied. We conclude that 
the essential physics of the velocity-dominated singularity can be captured by 
looking at homogeneous anisotropic spacetimes. 

Einstein's equations are contained in the lowest order /3-function equations of 
string theory (Lovelace 1986, Callan et.al. 1985, Fradkin and Tseytlin 1985, Sen 
1985); one does not expect the Schwarzschild solution to describe the geometry 
as one approaches the singularity, where higher order terms in the /3-function 
equations become important. Nonetheless it is worth repeating that there are 
exact conformal theories which contain spacelike singularities, including the 
SL{2,TZ) /U{1) black hole (Witten 1991), various four dimensional black hole 
backgrounds consisting of the tensor product of this theory with other conformal 
field theories (Giddings et.al. 1993, 1994, Johnson 1994, Lowe and Strominger 
1994), and a four dimensional cosmological model, with initial and final singular- 
ities, constructed as a c=4 WZW model (Nappi and Witten 1992). Furthermore, 
there are some indications that the SL{2,TZ)/U{1) black hole may be relevant 
for understanding the background near the singularity in higher-dimensional 
spacetimes; we will review these indications in the appendix. 

We will start with something rather simpler than these exact conformal 
field theories, looking instead for homogeneous solutions to the one-loop beta 
functions for the metric coupled to the dilaton (f> and matter with central charge 

2 

C = - (dcrit - D - A'^compact.) , (7) 

where dcrit is the critical dimension of the string theory, D is the number of 
macroscopic dimensions, and iVcompact. is the number of compactified dimen- 
sions. If we write the string frame metric as 

ds^ ^ -dt^ ' (8) 

i 

then the /3-function equations can be written as (Tseytlin and Vafa 1992) 

2 



D-i 2 / \ 

i=l V i=l I 

Afc - (^20 - ^ A,^ Afe 



2^-E^^-E^' = ^'"i- (9) 

i=l i=l 
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Here p is the energy density, and pk is the pressure in the direction of x'' : these 
terms come from adding a phenomenological matter action to the low energy 
effective action that can be derived from the usual /3-function equations. We 
may also include the effect of potential terms, such as that coming from a dilaton 
potential, in p and p. In general a will depend on how the matter couples to 
the dilaton. For a classical fundamental string source, the total action without 
the antisymmetric tensor field is (Dabholkar and Harvey 1989, Dabholkar et.al. 
1990) 

— 'S'spacctimc ~l~ 'S'worldshcct 

[X] , (10) 

where TZ is the spacetime curvature, TZ^^'' is the worldsheet curvature, and 
X{(t,t) is the location of the worldsheet. Here the fundamental string cou- 
ples as a classical source, and thus lacks the exp(— 20) = l/g^ij. factor. (We 
will discuss the justification for Equation (|l^) at the end of the following sec- 
tion.) Naturally, for multiple strings, one may sum the worldsheet action over 
the different sources, making the replacement 

5'worldshcct 

[X] > '5'worldshoct [^z] (H) 



in Equation (p^). Now if the stress tensor arises from a gas of classical strings, we 
expect that a = 2 due to the lack of a dilaton prefactor in the worldsheet action 
(Veneziano 1991, Tseytlin 1992). We will also get different values of a depending 
on whether the matter comes from higher loop effects and so on (Tseytlin 1992). 
In the particle production scenario we are proposing, the number density will be 
proportional to exp(— const./A), so that we will set a = 2 and let p = p{4>, A, A) 
and p = p{4'i It is worth noting that if we set the dilaton to be constant 

then Equation (|9|) has solutions only if 

D-l 

p=Y,Pk + 2c. (12) 

In particular, for c = a constant dilaton and a homogeneous metric requires a 
radiation-like equation of state. 

Many solutions to Equation (^ have been described by Tseytlin (1992): we 
would like to prepare for later discussions by sketching a few possible solutions. 
Define a new field ip — 20 — ^"^j^ A^ , and set a = 2; then Equation (|^) becomes 



4=1 



^2 ^ e'^p 
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n 
i=l 

(Tseytlin and Vafa, 1992). The first equation may be used to solve for (p, so 
that the second equation becomes 



D 

pe2^ + ^ A| - c (14) 

i=l 

With p,p — the solution will be a generalization of the Kasner solution with 

(pit) = ifio - ln(tc - t) 
Aj(t) = Xifl + qilii(tc - t) 

E'?' = 1 (15) 

1=1 

(Tseytlin 1992). Note that this corresponds to the + branch of the square 
root in Equation (|l4[). We wish to match this solution to the Schwarzschild 
black hole near the singularity, which as we have argued can be approximated 
by a Kasner metric. Thus, we would start with the solution to Equation (^5|) 
corresponding to (/> — 0. At this point we assume that particle production 
begins; what happens next depends on the details of the equation of state. In 
particular note that in order to turn around runaway expansion or contraction 
it is necessary that the pressure term become sufficiently large and negative 
for rapidly expanding spacetimes and sufficiently large and positive for rapidly 
contracting spacetimes. We will argue in the next section that string theory 
gives us exactly this sort of behavior. 

It is worth repeating here a point made by Tseytlin and Vafa (1992). In 
Einstein gravity we are used to any energy density at all contributing to the 
expansion of the universe, at least for a spatially flat cosmology where = 
(87rG/3)p. We are arguing that a negative pressure may in fact provide a drag 
on the expansion. The difference is the presence of the dilaton, which is an 
embarrassment for phenomenology but may be a blessing when dealing with 
singularities. Note also that runaway expansion or even a sufficient ammount 
of (positive) energy density will cause the dilaton to run to infinite coupling 
in finite time, making string loops especially important. On the other hand, 
sufficient contraction and insufficient energy density can cause the dilaton to 
run to weak coupling. 



(13) 



Afc = ^Pke"* ± Afc 



\ 
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3 String matter as a source for the field equa- 
tions 



Sanchez and Veneziano (1990) demonstrated that in a large class of inflationary 
spacetimes (with a constant dilaton), certain string modes became unstable and 
begin growing with the scale factor. They argued that this happens when the 
string size becomes larger than the horizon size of the expansion, so that different 
parts of the string cease to communicate with each other. They showed that 
this occured potentially with power-law expansion and de Sitter expansion, and 
inevitably with superinflationary expansion. Gasperini et.al. (1991a) studied 
the approximate behavior of large strings in the homogeneous, isotropic metric 

ds' = ^ (^-dv' + j:{dxr^ , (16) 

which describes power-law inflation for a > 1, de Sitter inflation for a = 1, and 
superinflation for < a < 1. They found that for sufficiently large strings, the 
conformal time rj was proportional to the worldshcct time t in conformal gauge, 
and that one could expand the solution to the classical equations of motion and 
Virasoro constraints in a power series expansion in the worldsheet time r and 
thus show that a consistent lowest-order solution to the equation of motion plus 
constraints is 

J? = X(a) + 0(t2) + ... 

V = T^f¥{a) + 0{T^) , (17) 

where < ct < 27r is the worldsheet spatial coordinate. Note that in this 

approximation, X' ^ X, and r) ^ r]' as r, rj — > (where the dot and prime 
refer to worldsheet time and space derivatives respectively), which is what one 
would expect from a largo, slow string. In a follow-up paper (Gasperini et.al. 
1991b), they showed that in rapidly contracting spacetimes one finds unstable 
behavior for which the worldsheet time is proportional to the comoving time t; 

the unstable behavior is characterized by X' <C X, and i ^ t'. (the former 
condition came from a term of the form t^~°'X2, where a < 1, which has a 
large time derivative as r — + 0.) 

Instead of repeating the calculations described above, we shall instead per- 
form a similar calculation to find the approximate behavior of very large strings 
in a temporarily expanding universe by using the string size, rather than the 
temporal distance to an infinite scale factor, as the expansion parameter. This 
will give us some intuition for the calculations we will perform in the following 
two sections. Consider an isotropically expanding or contracting universe with 
metric 

ds^ = C{r)) {^-dri^ + rfX^) . (18) 
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The classical equations of motion and the classical constraints in conformal 
gauge for a string in this spacetime are: 



-V + V - 



2C 







r]X - rj'X' 



C 



c 



f+rl'^-X -X^ 



--n-n' + X-X' 







(19) 



The first two equations are the equations of motion and the last two are the 
Virasoro constraints. We wish to expand these equations in a power scries in the 
size of the string. Let L be the dimensionless length, i. e. the full length of the 
string divided by igtr- Following Gasperini et.al. (1991a), we assume that the 
time dependence of the spatial coordinates is of lower order. We thus expand 
X in the series: 



X{t, a) = LXo{a) + i'""^n(r, a) , 



(20) 



where the Jf„ are all of order istr- We must be a little careful about the the 
worldsheet time interval we are interested in. As we will see, the lowest order 

solution to the constraints is rjo = tL\J X'q\ we will be interested below in 
universes which expand for a target space time interval of order where 
istrL/ p 1, so that the worldsheet time interval of interest will be of order 
1/L. In order to keep track of these factors explicitly let us define s = Lt, so 
that s will be of order 1 in this expansion. Then we may expand t] in the series 



(21) 



n=0 



In general we should note that we must expand C{r]) in a power series in 1/L as 
well, though this will not show up in the lowest order equations. These lowest 
order equations are 



iio + 



Xi 



P'2 



f,'o+X,-Xo 







cim) 



VoXi = 



jCi-X'o = 0, 



(22) 
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where the dot denotes a derivative with respect to s. One may solve these 
equations for arbitrary Xo{(j) with Xi = (so that we may absorb Xi into Xq 
and thus set Xi = 0) and 




The interested reader may check these equations to the next order and verify 
that solutions exist given these first order solutions. 

Note that in this case, just as in the case of the approach to the ry 0^ in 
the work of Gasperini et.al. (1991), the comoving string coordinates are frozen to 
lowest order (the time derivatives of X are of order l/L), so that the proper size 
of string will grow with the scale factor. We shall find that the same happens 
to the quantum wave functions below. 

If a large number of such unstable strings are produced as we will claim, we 
would like to know their equation of state, i.e. what the spacetime stress-energy 
tensor is for this collection of strings. For a single macroscopic string source, 
the action for the background fields coupled to this source is given by Equation 
(^0|). The stress tensor is then found by taking the derivative of this action with 
respect to the spacetime metric in the usual manner (Dabholkar and Harvey 
1989, Dabholkar et.al. 1990): 

T^,(x) = --^ / d^aVhh'"'dc.X^d,3X,e'''^S'' {X{t, a) - x)) (24) 
y G J 

Using this ansatz for the stress tensor, Gasperini et.al. (1991b) show that the 
approximate equation of state for unstable strings shrinking with the scale factor 
in a contracting universe is 

D-1 

which is the equation of state for photon- like particles; for unstable strings 
growing with the scale factor in an expanding universe (1991a) they find that 

D-1 

i=l 

This is the sort of behavior we desired in the previous section, in order to 
remove the singularity, although for general anisotropic spacctimes we expect 
the equation of state to be more complicated. 

The ansatz shown in Equation ( p^ for the coupling of a fundamental string 
source to the spacetime action still has the status of a conjecture, although 
there is evidence which supports it. Callan and Khuri (1991), Khuri (1993) and 
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Gauntlett et.al. (1994) have found that usmg this expression and calculating the 
low-energy scattering of two such strings using collective coordinate methods, 
one can reproduce the scattering amplitudes found from the low-energy limit 
of the Virasoro-Shapiro amphtudes. Tseytlin (1990) has attempted to derive 
Equation (|l^) by summing up long thin handles on the string worldsheet (these 
are essentially wormholes on the worldsheet). This is a very attractive and 
probably correct answer but it requires that one make a connection between the 
spatial size of the macroscopic string and the dominance of "thin handles" in 
the integral over the moduli space of high-genus Riemann surfaces. We do not 
understand how to connect spacetime scales to scales in these moduli spaces, and 
believe this to be an interesting problem. Still, the intuition that large strings 
should be treated as sources of spacetime energy-momentum, and that this effect 
should arise from the summation of string loops describing the interaction of 
the macroscopic string with its surroundings, is physically fairly sound; thus, if 
the effects of these strings are significant, we know that the dynamics of string 
theory causes certain classes of loop diagrams to become important, so that we 
must at least modify the /3-function equations appropriately. (Note that this 
is independent of another reason that loops might become important, namely 
that the dilaton may diverge near the singularity.) We shall comment on this 
point further in the conclusion. 

4 String field theory in curved spacetime 

4.1 How do we evolve the background forward? 

Naturally, we would prefer to understand the full quantum mechanical behavior 
of string field theory near the singularity, where field strengths and quantum 
fluctuations are large; we could then evolve the string field theory wave func- 
tional corresponding to the black hole background from a point where we trust 
our classical understanding of the geometry to the dangerous region near the 
classical singularity. Unfortunately this computation would require a mani- 
festly background independent quantum formulation of string field theory, not 
to mention a powerful computational control of that theory, both of which we 
lack. Even in the case of ordinary field theory, if we had a quantum theory 
of gravity, we would still have to contend with questions of interpreting the 
strongly fluctuating spacetime near the singularity. We will thus adopt the 
standard approach of computing the quantum stress tensor of matter in a fixed 
background, in order to use it as a source term for the /?-function equations of 
string theory. 

As we stated in the Introduction, this may be a reasonable approach for 
us. Since we are given as an adjustable parameter £pi/£str, we will perform our 
calculation in a regime where this ratio is small and we can trust the results of 
perturbative string theory; then, string physics becomes important well before 
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quantum fluctuations of spacetime become important. This is not necessarily 
the end of the story; if the background is composed of strings, then at the string 
scale there is no clean separation between the background and the string fields 
we are quantizing on top of it. 

In models with rapidly expanding directions there is another problem with 
keeping the background separate; if gravitons become unstable and grow with 
the scale factor, then the background gravitons become macroscopic strings and 
we must discuss the background as a soup of macroscopic strings. Perhaps 
there is some way of including this effect self-consistently to get some sort of 
modified /3-function equations which are valid in these regimes. Based on our 
discussion of macroscopic strings in the previous section, this would amount to 
the inclusion of a certain class of loop corrections to the /3-function equations. 
For the present, however, we will not attack this problem; rather, we will work 
in a background with a Hubble volume much larger than one in string units, 
and extrapolate the result to rjjubbic = -^string- If we can do nothing else, we can 
at least argue that the above issues must eventually be addressed. 

What we would like to calculate is the renormalized stress-energy tensor for 
string fields in a singular background. This is difficult in particle field theory, 
and should be at least as difficult in string field theory. As we have mentioned 
previously, in string theory we have at our disposal a few exact conformal field 
theories with spacelike singularities (Witten 1991, Giddings et.al. 1993,1994, 
Johnson 1994, Lowe and Strominger 1994, Nappi and Witten 1992). Since 
these examples correspond to gauged WZW models we might have some hope 
of calculating (T^'') exactly. However, several problems arise. The most natural 
calculation to perform is that of the one-loop expectation value 

Klcou^jx)) = ^^^^^^i J d^av'^^d^X'^d,X'^S{Xir,a) - x)) (27) 

in the first-quantized picture (this expression is derived by taking the derivative 
of the first-quantized one-loop partition function - which should correspond to 
the spacetime one-loop effective action for the background fields (Fradkin and 
Tseytlin 1985)). For this calculation we must know the expectation value for 
the general off-shell graviton vertex operator, and so we require an off-shell un- 
derstanding of the theory that we do not now possess. In addition, we must 
ensure that we are computing this expectation value (^^ so that it corresponds 
to the appropriate second-quantized vacuum, which means we must find some 
sort of prescription for dealing with the poles in the string propagator (recall 
that in particle field theory the vacuum is defined by the Feynman propagator - 
see Fulling's book (Fulling 1989) for a discussion). Finally, the quantization of 
strings propagating on non-compact WZW and gauged WZW models is not well 
understood. Given these problems, we will instead make a series of approxima- 
tions in order to make the calculation tractable. First, we will ignore concerns 
about the conformal invariance of the background, which will allow us to work 
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with simplified models of expanding spacetimes. Second, we will work only with 
the Lq constraint on the string wave functions, rather than the infinite number 
of equations arising from the full Virasoro constraints. Third, we will set the 
spacetime conformal time 77 to be a function solely of the worldsheet time r 
in conformal gauge (otherwise, solving for the string configuration at a given rj 
becomes extremely difficult). Finally, since we do not know how to compute the 
string energy-momentum tensor directly, we will compute the number of strings 
produced and from that calculation determine the stress-energy tensor. 

These approximations should not invalidate our results. First, we will be 
looking for semiclassical solutions to the string wave equation. As in the com- 
putation of the /^-functions using the background field method, when expanding 
the first-quantized path integral of the string around a classical solution, viola- 
tions of conformal invariance appear as logarithmic divergences in the quadratic 
fluctuation determinant. The effects we shall be interested in, namely the ex- 
ponential dependence of the production amplitude on the energy, come from 
the classical term exp{{i/a')Sci} of the wave function. Second, although we 
use only the Lq constraints rather the full set of positive-mode Virasoro con- 
straints, imposing the latter should simply have the effect of reducing the degrees 
of freedom of the strings that we produce. A useful analogy is the computation 
of Hawking radiation in QED; if we work in the Lorentz gauge and forget about 
the Gupta-Bleuler quantization procedure, we get the same basic physics out of 
the calculation, but we find that the calculated Hawking flux is too large by a 
factor oi D/{D — 2), due to the effect of the unphysical timelike and longitudinal 
modes of the photon. A similar story should describe our calculation as well. 
Third, although setting i] = ri{T) is not consistent with the full Virasoro condi- 
tions, as we have seen, an examination of Equation (|l^) shows that this ansatz 
is consistent with the Lq constraints. Finally, although generally the concept of 
single particles or single strings is ambiguous in curved spacetime backgrounds, 
we will work with a spacetime in which string number does make sense, namely 
a temporarily expanding model which becomes flat Minkowski space in the far 
past and in the far future. Here we may sensibly calculate the number of strings 
using standard techniques and then simply multiply the number of strings in 
a given quantum mechanical state by the spacetime energy-momentum tensor 
of a string in that state, using the results described in the previous section, or 
some modification of them. Our methods will not be directly extendible to a 
calculation of the quantum stress tensor in a spacetime in which the future is a 
singularity rather than Minkowski space, but it will allow us to make some sen- 
sible statements about quantum string field theory in curved spacetime which 
should apply to the region we are interested in. 

4.2 Imposing the Lq constraints 

We wish to calculate the number of different strings produced just as we would 
calculate the number of particles produced: by finding the solutions to the first- 
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quantized relativistic wave equation which are purely positive frequency in the 
far past and looking for the coefficients of negative-frequency components in the 
far future. For string theory, the (linearized) wave equation is just the Virasoro 
condition: 

r++ + r__ = G^^(x(a)) (xA^x^ + x'^x''^) =0 

T++ - T__ = G^^X'^X'' = . (28) 

We wish to turn the Virasoro operators into functional differential operators 
acting on wave functionals '9{X{(7)). In curved space, there are subtleties even 
in the case of particle wave functions due to the necessity of defining the Hilbert 
space norms and the path integral measure in a covariant fashion (DeWitt 1957); 
we will assume the same story applies when we derive the string wave operator 
from the first-quantized path integral. The functional differential operator form 
of (M) should look like 



T++ + T__ 



G^,{X{a))x'^'X''' +CH 
T++-T__=X'^- , (29) 

where TZ is the spacetime Ricci scalar and ^ depends on whether the string is 
conformally, or otherwise, coupled (in other words, with what coefficient we 
multiply a term like TZ{X{a)) in the first-quantized string Lagrangian.) We 
prefer to remain agnostic as to how the string is coupled; at any rate, in the 
approximations we shall look at the curvature term will be of lower order and 
so we will ignore it. 



There are several issues we must deal with before using Equation ( |29[ ) as 
a wave operator. First, it must be normal ordered. Secondly, it represents 
an infinite number of differential equations, one for each Fourier mode of each 
operator. In general, these Fourier modes are quite complicated; even if G 
depends only on the conformal time 77, 77 is expandable in Fourier modes and the 
metric is generally a nontrivial function. Finally, it is difficult to understand and 
interpret the wave functional in spacetime since the spacetime time coordinate 
of the string depends on the worldsheet spatial coordinate, as it does in Equation 
(23). With this in mind we will use only the ct- independent modes of Equation 
(29) (in other words, the Lq and Lq constraints), and we will set rj — rj^r). 

If we write the Fourier expansion of our string modes as 



00 



X^{T,a)^ > X„(r)e™'^ (30) 
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and use the reality constraints, we find that Xn = ^In- Let 

x';, = xl^ + ixl^ , (31) 



so that X_n = X^ — iX^. We then find that the Lq part of ( P9[ ) is, for the 
homogeneous spacetimes we are interested in, 

rf-1 



k=l 71=1 



d-1 / o \ 2 



k=l 
d— 1 C30 



fc=l n=l 



^2,„T7^ 1 !> {^n}]=0, (32) 

(where we take the metric to be diagonal). Here ^^jl^n = d/dXn^^. Imposing 
the final constraint in Equation (|3|) becomes much simpler if we set 

X„ = R^e'^- , (33) 

as Equation ( [32| ) becomes: 



r*^*^ ( ^ ^ u'' ^ 1/9^^ 

d—1 OD \ 

YY''' (^^)'-e7^ U[r;,{X„}]=0 

k=ln^l ) 
d—1 oo ^ 

EE'^'M^*['?'{^»''^"}]=0. (34) 

fe=ln=l 

Note in particular that these constraints commute with d/d(f)^^ for all n,fc, so 
that we may let 

* [{R, cj,}] = exp{z Y ^fc,n</'n}* [{^}] , (35) 
and thus make the substitution 

d 

TTTT ■ (36) 
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It will be useful at this point to write down the stationary states of ( p4[ ) 
in flat space. We will add a term £q'^ to this equation in order to explicitly 
subtract off the zero-point energy: this will include the usual tachyon shift 
(— rrig) if we are working with the bosonic string. In flat space the equation 
clearly separates into an infinite number of 2-dimensional harmonic oscillators, 
one for each {k,n), with the single constraint that 

d—l oo 
fc=l n=l 

We also get a Klein-Gordon equation for the Xq wave functions. The final wave 
function may then be written in terms of associated Laguerre polynomials. If 
we leave the scale factor of the metric in explicitly, we find that 

M/r , 1 = -iE-n+iq-Xo „ 



He 

k.rik 



(38) 



Here 



rrik.n + i ) - f o' : (39) 



and the final line in ( p^ ) comes from (|3j). Note that acts as the frequency 
of the harmonic oscillator, and the only difference between this system and a 
product of nonrelativistic harmonic oscillators is that E gets replaced by E^. 
As usual, Sq is defined such that 

=q^ +Y1 4"™fc,« - ml . (40) 

n,k 

4.3 Approximations for cosmological models 

Let 

d-l 

ds^ = -dt^ + ^ ak{t)dxl . (41) 
fc=i 



Putting this metric into Equation (34), we find that the scale factors a{t) act 



as time-dependent masses for the different 2D harmonic oscillators. In this 
case. Equation is not separable, and no interesting scale factors yielding 
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exact solutions are known to us. (This is in contrast to the non-relativistic 
harmonic osciUator with a time-dependent mass and/or frequency, where one 
can reduce the quantum problem to the classical problem by introducing a 
Gaussian wave packet - essentially a coherent-state wavepacket ~ with time- 
dependent coefhcients (Peremolov and Popov 1969). There is no solution in 
the relativistic case as there are no analogs of coherent states: even with a 
time-independent mass and frequency, the relativistic dispersion relation causes 
the wave packet to spread.) We will instead examine Equation (|3^ ) in certain 
limiting cases. 



4.3.1 The adiabatic limit 

The simplest limit to study is the limit of an adiabatically changing universe, 
where ak = ak{pt), l/p <C istr (note we have so far set £str = !)■ Let us work to 
lowest adiabatic order; since the curvature scalar is second order in derivatives 
of the metric, and thus 0{p'^), we will ignore this term for the present. To this 
order, the oscillator number should be an adiabatic invariant, so we use the 
following variant of the flat space energy eigenfunctions: 



n 

k.rih. 



(42) 



If we now expand the total wave function in this basis using time dependent 
coefhcients, 

E q™.,.,'...}W*{™.,„A,„} ' (43) 

{mk,„,lk,n} 

then applying Equation (^J) to this wave function will give us a coupled set 
of differential equations for the coefficients Cj^j. ^ j^, „}(i), subject to the initial 
conditions appropriate to the question one is asking. 

We would like to compute the probability of particle production. As usual, 
we first need to define the positive frequency modes of the string field in the far 
past and future. The standard definition is the definition of adiabatic positive 
frequency modes (see Birrell and Davies (1982) for an explanation and refer- 
ences): we ask that the modes we define as "positive frequency" be those that 
match exactly onto positive frequency WKB solutions in the infinite past, 

*S.„,.,„,;^»P{-/_J«(0} , (44) 

where for our spacetime 



^('') = \ / E ^ + E + . (45) 
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Ideally, we would find the appropriate exact wave function that matches onto 
this WKB solution and look for the negative frequency components in the usual 
manner: in lieu of this, we would like a WKB approximation for the wave 
function at all times. Two problems arise, however. The first problem is that 
the negative frequency components will not be present at any order in the 
perturbation series one would get from plugging ( ^3|) into (^J). The reason is 
identical to the reason one does not find above-barrier reflection in scmiclassical 
non-relativistic quantum mechanics to any perturbative order in h. The second 
problem is that the terms in the equations for the coefficients Cm,i{t) in Equation 
( ^3|) which are higher order in p come from applying the time derivatives to vj/^'^-^ 
which depend implicitly on time through au ■ The reader may easily verify that 
these terms are also multiplied by factors of mk,n and higher powers of 771, 
because of the recursion relations for Laguerre polynomials. Thus if mkjp 3> 1, 
the adiabatic approximation breaks down. We will treat this case below. 

In the case that the adiabatic approximation holds, there is still a trick for 
extracting the negative frequency part of the wave function in the far future. 
Audretsch (1979) has pointed out, for scalar particles in curved spacetime, that 
the particle production calculation is almost identical to the calculation of above- 
barrier reflection in nonrelativistic quantum mechanics, since one can generally 
(for homogeneous spacetimes) separate the time and space parts of the wave 
equation and write the former as 

^/ + c.(t)V = 0. (46) 

The difference is that in above-barrier reflection, we want the incident compo- 
nent of the wave function to be normalized to unity and the transmitted wave 
function to be purely outgoing. Then, if the amplitude of the transmitted wave 
is T and the amplitude of the reflected wave is R (assuming that only one value 
of jwt^ool contributes to the final wave function, though the discussion may be 
easily generalized) probability conservation requires that 

|r|2 + |i?|2 = 1 . (47) 

On the other hand, in the case of non-static spacetimes, we wish the wave 
function in the far past to be purely positive frequency. If this incoming wave 
function is normalized to unity, and the coefficients of the positive and negative 
frequency components of the outgoing wave function are denoted a and P (again 
we assume that only one magnitude of the frequency contributes to the final 
wave function), probability conservation requires that 

1 + = lap . (48) 

If we write ui = \/2{E — V{t)) the analogy is obvious. In general, oj will vanish 
in the upper half complex t-plane at some point (or set of points) to, and also 
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at Iq. The square of the reflection coefficient in the WKB approximation has 
been derived by Pokrovskii and Khalatnikov (1961): 

= exp |-2 Im^^ dt cj(<)| . (49) 

In our case, if we substitute E for oj, becomes the square of the Bogolubov 
coefficient multiplying the negative frequency WKB solution in the far future 
(Audretsch 1979). 




Re t 



t * 

Figure 1 

Complex time paths, to and are the classical turning points of the system. 
Curve C denotes a path which moves around the square root branch point, 
while staying far enough away for the adiabatic or semiclassical approximation 
to hold. For the purposes of computing the production amplitude the integral 
J dtE may be performed along curve C. 

Before we continue, let us describe a more heuristic derivation of Equation 
( p9[ ) (see also Landau and Lifshitz (1977) and Migdal (1977), from which this 
discussion was taken). As we have seen, the time dependence of our system 
has the form of Equation (^). The expression for E in equation (|4|) contains 
square root branch cuts in the complex t plane (see Figure 1). In order to detect 
the part of the wave function with negative energy, we should go around the 
end of the branch cut in E. If we stay far enough away from this endpoint, 
then the adiabatic approximation will remain valid. (The approximation will 
break down near the turning point, where E = 0, since one requires for the 
approximation that the time variation of the wave function be more rapid than 
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the time variation of the background.) Along the contour, labeled C in Figure 
1, the integral / dtE will pick up an imaginary component, which will give us 
the real part of the Bogolubov coefficient f3. This real part will not change if we 
move the points al,a2, where the contour C leaves the real axis, adjacent each 
side of the branch cut; so long as C does not enclose any singularities other 
than the branch cut, we may then shrink the contour down to the curve C 
which surrounds the branch cut as shown. A little manipulation of the integral 
should convince the reader that the imaginary part of the integral in Equation 
( ^9[ ) is then the correct result for the real part of the phase picked up when 
traveling on contour C. Since exjp{—E/p} will be small in the adiabatic limit, 
the correction for the correct normalization of the wave function will be even 
smaller, and we may identify the expression shown in equation ( ^9|) with 
This discussion can be easily generalized to systems with several turning points 
in the complex plane (Pokrovskii and Khalatnikov 1961, Landau and Lifshitz 
1977, Migdal 1977, Audretsch 1979). 
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Figure 2 

Pair production of strings. Here we depict the configuration of the classical 
string as it passes around the branch point in the complex t plane. One may 
interpret this as a pair production process beginning at to; the string which 
reaches the real axis on the left of the branch cut is positive frequency and the 
string which reaches the right side of the branch cut is negative frequency. 

One can envision the pair production process via the classical trajectory of 
the string on C . Since the energy is small near the turning point, the size of the 
string will be small as well. The classical trajectory of the string describes two 
zero-size, zero-energy strings being produced at to, and growing as they travel 
in imaginary time along opposite sides of the branch cut (see Figure 2). 
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4.3.2 The semiclassical limit 



In the limit that the occupation numbers mk,n are large, the adiabatic approxi- 
mation breaks down and the semiclassical approximation is appropriate. Let us 
write the wave function as 



'S,(t,{^„}) 



(50) 



In the semiclassical limit, where the frequencies (excitation numbers and center 
of mass momentum) are large, S will solve the Hamilton-Jacobi equation 



dS_ 
'dt 



E 



k\2 



O-k 



^ akit) 



k.7i 



/ dS 







dS \ 






d4>U 



- = 



(51) 



where S defines the momentum as a function of position via p„ = dS/dqm p is 
a multivalued function defined by a curve C in phase space (see Figure 3) and 
the sum in (js^) is the sum over the branches of the function p{q) defined by C. 
If we foliate phase space with curves C ~ C{'P), then we may write Dj as the 
square root of the Van Vleck determinant 



Dj =^ constant x 



dqidVk 



(52) 



(This geometric picture is taken from Berry and Balasz (1979) and is due to 
Maslov (1965). See Berry (1981) for further review and references.) In general 
the condition for the semiclassical approximation to hold is that: 



dql 



dS 
dqn 
dS_ 
'dt 



(53) 
(54) 



This is the usual criteria for the validity of the WKB approximation. 



24 



p 




q 

Figure 3 

Curves p{q) in phase space. The n- dimensional surface Cn{V) defines a 
multivalued function pk{{qi}) at time rj in a 2n- dimensional phase space. V 
labels a family of n- dimensional curves which foliate phase space. C^i(7^i) is 
found by evolving each point on C^(7^i) from time rj to time rj' via Hamilton's 
equations. 



Note that in general the Hamilton- Jacobi equation for a relativistic system 
should look like 



as 
a7 



n[t,x,E 



dS 
'dt 



,P - 



dx ' 







(55) 



if we are to treat space and time on an equal footing, r should be an arbitrary 
parameter describing the world line of this system, i.e. the r dependence of the 
system should be trivial. Thus, we arrive at Equation (51), where 



-H{t,E,Rtp'k,n,€,Pin) 



E 

k 



k\2 



Ml 



E 

k,7i 



(Pr.u) 
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ak{t) {RiYakit) 



ak{t){K) 
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+C7^ - mj; = . (56) 

The simplest way forward is to simply solve for the constraints, by letting E, the 
generator of t translations, become the Hamiltonian and using Equation ( ^ ) to 
write 



(57) 



This also comes from putting the term (dS/dt)'^ on the right hand side of Equa- 
tion (^) and taking the square root of the equation. 
Now we must solve the Hamilton- Jacobi equation 

f + i^^O. (58) 

We describe the solution as follows (this discussion is patterned after that in 
Berry and Balazs (1979), though our representation for S is different). At some 
initial time, say t = 0, we start with a curve Cq in phase space which defines 
the desired initial wave function via (|50|). This curve will evolve forward in 
time, as one can see by taking each point on the curve as an initial condition 
for Hamilton's equations, to a curve Ct at time t. This will define a family of 
functions p{q,t). Assuming no new caustics in p{q,t) develop, the solution to 
the Hamilton- Jacobi equation can be written, to within a constant, as: 

5= / Pn{q,to)dqn - I E{qn,Pn{qn,t'),t')dt' . (59) 

Here the q{n,o) are arbitrary constants. A straightforward application of Hamil- 
ton's equations confirms that 

and that S as defined in Equation ( |59| ) satisfies the Hamilton- Jacobi equations. 

One may compute string production in the same fashion as one would in 
the adiabatic approximation: here one would search for turning points in the 
complex time plane for p„ and E. To see this procedure in an explicit model 
we turn to the next section. 



5 An isotropically expanding model 

We would now like to examine a model from which we may extract a few results. 
For the reasons laid out in the beginning of the previous section, we would like to 
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examine a model in which the far past and far future correspond to Minkowski 
space and there is a finite period of expansion near 77 = 0. At present the one 
case we have found to be tractable is the case of isotropic expansion. In this 
case we will work in conformal time: 

ds^ = drj) l^-dT]^ + j . (61) 

A simple scale factor with the desired properties is C{ri) = A + B tanh pr], where 
A — B — 1 and A+B — fl, the final scale factor. The conformally coupled Klein- 
Gordon equation in a two-dimensional spacetime with this scale factor has been 
solved exactly and examined by Bernard and Duncan (1977); the time part of 
the equation is identical to the Schrodinger equation in a special case of the 
Eckhardt-Sauter potential (Eckhardt 1930, Sauter 1932). It will be useful as 
a warmup exercise to compute the production probability for particles in this 
spacetime using the WKB approximation, since we have an exact solution to 
compare it to. After some work we will find the semiclassical result for string 
production to be nearly identical. 

5.1 Review of the Klein-Gordon equation 

For a particle with mass m, we are looking for solutions which look in the far 
past and far future like: 

t — * -00 : (l){x, t) - e-''^-*+*«^ 
t^oo: 4>{x, t) - a,e-*'"+*+*«" -|- /3,e''"+*+''" , (62) 

where 

ujl=q^ + {A± B)m^ . (63) 

Bernard and Duncan (1977) found exact solutions to the Klein-Gordon equation 
(in the form of hypergeometric functions) which match to (|62[). The exact 
expression for the production probability is 

sinh (vrcj- / p) sinh {ttuj^ / p)) 

In the limit of small p (or large energies and energy differences) , this expression 
becomes 

= e-2---/p . (65) 



Using the WKB approximation, wc find that Equation (^9|) becomes 

|/3fcp =exp|-2Im^' y/k^ + {A + B tanh pt)m'^^ , (66) 



27 



where 



in 1 ^ / cj- 
2p 2p 



Lu^ \ In z 



(67) 



(68) 



(69) 



Plugging this into Equation (|66|) we find that this approximation reproduces 
Equation (BSh. 



If we make the substitution 

r/ = In 

then the integral in Equation (|6^) becomes 



LO- f dz Z + 1 ITTUJ- 



5.2 Strings in an expanding universe 

For states for which the adiabatic approximation holds, it is easy to see that if 
we ignore (^7?, — mg) in the equation, we may use the approximations described 
to arrive at Equation (|65|). However, states for which this approximation holds 
are not generic at sufficiently high mass levels. String theory contains an infinite 
number of states with a multiplicity growing as the mass. For a given level M = 
(mass)^ of string excitation, states corresponding to high Fourier components 
K of the string excited to low oscillator levels rriK, such that rriKP 1, can 
still be treated with the adiabatic method. But in general the typical size of 
the string will be large, reflecting a large number of Fourier components excited 
to a high oscillator level. We will find that while some distributions of the 
energy for a given M are intractable even in the semiclassical limit, i.e. those 
with a few Fourier components highly excited, a configuration with N ^ 1 
Fourier modes excited may be treated in a 1/N expansion. The latter type 
of configuration should be fairly generic for sufficiently large mass; we expect 
that most configurations at a given total oscillator level n will have their energy 
partitioned among many oscillator modes. 

In conformal coordinates, the scale factor appears as a time-dependent fre- 
quency rather than a time-dependent mass; thus, the energy (the generator of 
r] translations) is 
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Hamilton's equations are: 
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(71) 

(There is also the remaimng constraint ^ nlk.n — 0.) These equations may 
be converted into equations of motion for R and (j): 



Ri, 



E' 



(72) 



Note that without the E/E terms, Equation ( [72| ) looks like the equations for a 
harmonic oscillator with frequency 



nC{r]) 
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If a single oscillator is excited, Equation ( [7^ becomes quite difficult to solve 
in the limit Ep/u) ^ 1. The E/E term is of lower adiabatic order, but the 
reduced frequency is smaller still. The period of oscillation in spacetime time 
is much longer than the time scale /5~^. Whether loR/E is of order 1, of order 
p, or of order and thus whether the velocity- or position-dependent force 
dominates, will depend on where in the potential the oscillator term is. At 
present, we have no ideas to offer for a solution in this case. 

However, as we have stated above, most of the states at a given large mass 
will be in oscillators with many Fourier modes excited - call this number N . 
We will attempt to solve for the classical equations to lowest order in This 
is similar to the expansion in the size of the string we presented two sections 
ago; the string in such a state will be large and slow. We will assume that for 
a given oscillator mode (fc,n) with energy i?fe^„, E^^nP/^ 2> 1 and N ^ 1; thus 
E > Ek^n for any (fc, n). 

Using Hamilton's equations we can see that 



E 
E 
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The expression multiplying C/C is 0(1) both in 1/N and in p. The position- 
dependent terms in ( [7^ ) are of order 1/N. Even though Np ^ 1, the oscillator 
terms will dominate for most of the life of the universe, where the tanh(p77) term 
in C is exponentially suppressed and so not changing much. More specifically. 
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note that 

C Bp sech^pr; 



C A + B tanh p?] 

- Bpe"2p|r,| foj. ^ — > ±00 . (75) 

Thus, noting that Hamihon's equations tell us that R ^ 0{1/^/N), so long as 

\^\>l-\n{\B\p^}=r^f , (76) 

the position dependent forces dominate and the system behaves as a collection 
of simple harmonic oscillators. Since 

i?~^o(VA^) , (77) 

for I77I > r]f the change in E (or in co) becomes very small exponentially rapidly, 
and the system behaves as set of oscillators with frequency 

_ nC{±oo) 

Note that the factor of means that the potential is quite shallow and in 
general each Fourier mode will have a spatial extent of order VTV : the string 
will be quite large, as promised. For \ri\ < rjf, the velocity dependent term 
dominates. Keeping only these terms and the second derivative in ([72|), we find 
that 



a'' 



E 



k = (79) 

where the a are integration constants and will be fixed by matching this solution 
to the I ?7 1 > 77/ region. Thus the velocities are 0(l/\/]V): we shall call this region 
the "frozen region" (thus the suffix in 77/). The changes in these coordinates 
are: 

5{R, cjy) ~ 2rjf{R, ^) ~ In {\B\pVn} . (80) 



pVN 

As long as p\/N ^ 1, the change in the coordinate positions is small; therefore, 
we will approximate the string in this region as being fixed in conioving coordi- 
nates, so that its proper size grows with the scale factor, as we saw two sections 
ago. 

Physically, the picture that we have is that strings which are larger than 
a horizon radius during expansion become "frozen" and expand with the scale 
factor without oscillating, as different parts of the string are no longer in causal 
contact with each other. The entire history of the string can be divided into 
three regions: 
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• —oo<r]<r]f. Strings oscillate freely and energy is conserved. Here we 
may expand the string wave function in stationary states of the initial 
Hamiltonian. 

• —rjf < rj < rjf (the "frozen" region.) Classically, strings larger than the 
horizon size do not oscillate coherently since different regions of the string 
are causally disconnected from each other. Instead, the string grows with 

the scale factor. Eqnivalcntly, in conformal coordinates the string wave 
function is frozen as the harmonic oscillator potentials decrease in width. 

• r]f < 7] < CO. The wave function at 77 = rjf is no longer a stationary 
state of the final Hamiltonian: we may, within our approximation, de- 
compose it into eigenstates of the final Hamiltonian, just as in the sudden 

approximation. 

In the language of Berry and Balazs (1979), if we are looking for a wave function 
corresponding to an energy eigenstate with energy E- in the far past, we wish 
to start with a curve C_oo m phase space corresponding to the torus defined by 
the equation E = E- in phase space. In this region the energy of each oscillator, 
i?fe_„, is a constant of motion, so that 

p>k^^ = ±^{E^r-u;UR'nr-{^^)' ■ (81) 

The branch of the square root changes when we reach a caustic in p{q) and, 
as stated, we must sum exp{iS} over the branches of p. In our approximation, 
when T] < —rif the energy is effectively time-independent and the resulting semi- 
classical wave function is just the semiclassical approximation to the product of 
harmonic oscillator wave functions. In the "frozen" region, since R, (p, and pn 
are 0{l/y/N), to lowest order in our 1/y/N expansion the curve Cjj will remain 
frozen. Within our approximation we may write the classical action as: 




The curvature term in the expression for E is 0(1 /TV) smaller than the other 
terms, so we shall ignore it. As we leave the frozen region, we may calculate 
the negative frequency component as before (the reader may easily check that 
the turning point is well within the region that the "frozen approximation" is 
valid), and we find once again that the Bogolubov coefficient is: 

|/3{„}|2 = e-2-^-/^ (83) 
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which is essentially identical to the expression given in ( |65| ) . This is independent 
of the position, despite the fact that / dr]E, and the turning point in the complex 
77 plane, depends on R. For more general scale factors the imaginary part of the 
integral around the turning point of E will depend on R. In addition the time 
integral in the "frozen" regime has some position dependence, so that the wave 
function begins at 77 = — 77^ as a stationary state and in the "frozen" region 
begins to develop additional structure. Nonetheless, within the approximation 
we are working in, the position dependence of Im J drjE and the perturbations of 
the initial wavefunction which develop as we move through the "frozen" region 
are subleading. One can see this by noting that 

dE {Cjrj)^ - 1) Rj 

ai4 " E ' ^^^^ 

thus the change in the energy arising from the change in any given oscillator 
coordinate will be 

E'iR^O{\l^). (85) 

For computing the integral Im J dr\E^ summing the squares of this error gives 
us a total error of 0(1), compared with a lowest order answer E_ / p which is 
0{N). In examining the perturbations to the stationary state wavefunction 
that accrue during the evolution through the "frozen" region, we can see that 
for a change in a single oscillator coordinate. 



■,jd,E^2,,SE 



(86) 



while the integral / dR'^p^ is of 0(1) in the expansion, so that corrections 
to our description of the wavefunction as being frozen in this regime are indeed 
subleading. 

After the string leaves the frozen region, the curve C^j, will generally intersect 
some band of energy shells of the Hamiltonian describing the evolution beyond 
77 > rjf, as shown in Figure 4. The curve Cri>nf becomes time dependent and 
starts to rotate and stretch with a periodicity given by that of the final oscillator 
Hamiltonian. 

We still need to check that our semiclassical approximation is valid, by check- 
ing Equations (|5^) and (p^). Hamilton- Jacobi theory tells us that dS/dq = p{q). 
Thus, the condition for Equation ( p3| ) to hold is that: 

« {vtf ■ (87) 

This is just the usual condition for the semiclassical approximation to hold for 
the harmonic oscillator wavefunction (though in the relativistic case this means 
uj'^q <C rather than uj'^q ^ p^)] it will work for large energies so long as we 
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are sufficiently far from the turning points. The energy is large and fixed for 
77 < ?// in our approximation; at rj = ijf the curve Crif intersects a finite band 
of energy shells of the final Hamiltonian, all with energy larger than the initial 
energy if the scale factor expands (and smaller if it contracts, though so long 
as the change in cu is not too drastic the final range of energies should still be 
large). For states with oscillator modes of a given energy E^^n, Pr is given by 
Equation (|8^), and thus Equation (|5^) holds so long as E^^n is large. In our 
approximation this is the case; if E^- n was small we could treat that oscillator 
via the adiabatic expansion. To check Equation ( |5^ ) , note that 

Equation (|^ holds if 

which is certainly the case. One may also work with the path integral of the 
wave function and check that the quadratic fluctuation operator, given by second 
variation of the action along the classical trajectory, contains terms of order E, 
so that the saddle point integration around the classical trajectories gives a good 
approximation. Recall also that it is the divergences in the determinant of this 
quadratic fluctuation operator which contribute to the /3-function of the 2-d a 
model; we would thus like to reiterate that although our background does not 
lead to a vanishing /3-function, the effects of this lack of conformal invariance 
will appear in the determinant of the quadratic fluctuation operator, while the 
physics we are interested in resides in the exponential part of the wave function, 
ex.p{iSci.}- 

For r] > rjf, the oscillator frequency has changed and the system now has 
frequency uj+. It is no longer in a state of definite energy, however. As shown 
in Figure 4, the curve C,,^, which is still in our approximation identical to C_oo, 
intersects many tori defined hy E = constant when u> = w+ and so Cri^nj will 
be time dependent. However, as we know that the wave function is stationary 
for rj < rjf within our approximation, it is easier to write the spatial part of 
the wave function as the product of Laguerre polynomials that our semiclassical 
wave function approximates; we may then simply decompose this wave function 
into energy eigenstates of the Hamiltonian for rj>rjf. The reader has doubtless 
noted that this is essentially the sudden approximation for the spatial part of 
the wave functions. We may argue for this approximation more directly using 
the Hartree approximation for the wave functions: that is, let 

^ = l[i^kARn,^tv) , (90) 

k,n 

and assume that correlations in the wave function are of higher order in 
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We may write the wave equation applied to Equation ( |90| ) 

J2nff''\Rtv)\^ = o (91) 



k.n 



as 



k,n I {l,m)^{k,n) {l,m)^{k,n) {j,p)^{k,n),{l,m) 

(i,m)^(fc,n) J 

The resulting equations of motion are 

9^ . . d 



where 



Q^f A^^fco,no(^)^ + Wfeo.no (fi, ^7)1 V'fco.no = (93) 
NKkAv) ^ I ^lm^^l,m ■ (94) 

Let us rescale rj N-q. The first term is now 0{\/N'^) and we shall drop it. 
The time dependences of K and of the frequencies in Ti have the form TV 77, so we 
have a non-relativistic harmonic oscillator which varies rapidly compared to the 
other time scales of the problem. In particular, C{Nrj)'^ — A + B ta,nh{N prj) is 
effectively a step function when Np 3> 1, so the sudden approximation is valid. 
This argument is admittedly very rough but seems to capture the right physics. 
A string with many oscillators highly excited will be heavy and slow; therefore, 
it will be nonrelativistic and any change in the spacetimc, if moderate, will 
occur on a time scale much smaller than the response time of the string. We 
should note that we still need the relativistic analysis above to find the string 
production amplitudes, as there is no string production in the nonrelativistic 
approximation. 
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Figure 4 

Evolution of one of the Fourier modes of the string. For rj < rij pk{qk) is 
defined by the solid curve C—QQ, which is a curve of constant energy in 
phase space. For ri < rjf the curves of constant energy are the dashed ellipses. 
The dotted curve Cn>nf denotes the function p{q) found by evolving C_oo from 
r]f to some r] using the string Hamiltonian. 

We have thus shown that in this model, the number of strings per unit volume 
produced in a given state with energy is proportional to cxp{—2nE^/p}. At 
sufficiently large oscillator level N, if the center of mass momentum is of order 
1, -B- K, m K, y/N/a' where m is the mass. It is well known that the number 
of string states of mass m, for m 1/a', can be approximated by the formula 

D{m) » m-''e^^" , (95) 

where b and c are pure positive numbers of order 10 and depend on the specific 
string theory one is using. So far we have calculated the production amplitude 
for an arbitrary state of a given mass regardless of whether or not it was a 
physical state. The density of states should automatically count physical states, 
so if we use it to count the total number of strings produced we should get the 
correct answer. Thus the total energy density produced during the expansion 
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is, if we assume that our production formula holds for most string states at a 
given mass, 



unitvolume 



/•oo 

/ dm mD{m)e-^''"'^P . (96) 
Jq 



If we allow ourselves the luxury of extrapolating our result to expansion rates 
of order the inverse string scale Vc/, we find that when 

P > ^ (97) 



the integral in ( Pq ) blows up in a manner reminiscent of the Hagedorn transition. 
This is the main result of this paper. 



5.3 Interpretation 

The fact that the energy produced at a critical value of the expansion rate 
diverges indicates to us that computing string production on a fixed background 
leads to a diverging rate of string production. We believe this divergence to be 
a sign that we have either reached some sort of phase transition in the theory, or 
we have neglected the effects of the backreaction of the produced string matter. 
We may speculate in the latter case that since the bulk of the energy is produced 
in very large strings as p approaches 1/Va', the backreaction should act to slow 
the expansion, as we have previously claimed should occur with a sufficient 
density of strings which grow with the scale factor. Since this problem will 
arise whenever the expansion rate is sufficiently fast, we argue that the effect 
of backreaction will be to sufficiently slow the expansion so that any curvature 
singularity, if reached, is reached at infinite proper time and that all geodesies 
become completeable. 



6 Speculations about anisotropically expanding 
models 

The arguments in the previous section are far from a rigorous solution to the 
problem of generic singularities, which as we have argued are very anisotropic 
and furthermore have no "out" region with an unambiguous definition of particle 
number. As for the latter, we hope that the general features of our answer in 
terms of the more covariant concept of the energy-momentum tensor will carry 
through to a calculation of this quantity in a background with diverging rates 
of expansion and contraction, even though the details of the calculation will 
be very different; in particular the essential singularity in the calculated rate 
of energy production seems physically well-motivated and general, being based 
on the fight between an exponential suppression of string production and the 
exponential growth of the density of states. 
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If we wish to understand the further evolution of the spacetime, we need 
to understand how the anisotropic expansion and contraction which exists near 
generic singularities affects the configuration of produced strings, and we need 
to understand the equation of state in these anisotropic spacetimes. For exam- 
ple, one might expect that strings are preferentially produced with excitations 
in the most rapidly expanding and contracting directions, and that most of the 
strings produced will in the expanding case be large strings growing as the scale 
factor and in the contracting case be small strings with large momentum in 
the oscillator modes. The wave functions we have constructed are stationary 
states with no definite shape. To make our scenario plausible we would need to 
work with some set of coherent states which describe strings in definite classical 
configurations; since the highly massive strings are generally nonrelativistic the 
width of the string wave function around the peak configuration will not spread 
too quickly. We have not been able to detect any preference in the produc- 
tion of any particular configuration at our level of approximation, even in the 
isotropically expanding case. 

As for the question of which string field modes are preferentially produced, 
if we examine strings which are polarized in a direction or set of directions with 
a given expansion rate, the number of strings produced will increase with the 
expansion rate in that direction. On the other hand, strings with no polariza- 
tion in the expanding and contracting directions (if some spatial directions have 
a constant scale factor) should not be produced at all. There should be some 
extrapolation between these extremes, but our approximations seem to be rela- 
tively insensitive to the distribution of modes among the various polarizations. 
Clearly there is a considerable amount of work to do. 

One natural scenario for the evolution of the spacetime is the following: as we 
approach the singularity, where it is approximated by the Kasner metric, string 
production becomes important, and we produce unstable strings preferentially 
polarized in the directions which are changing the most rapidly. The equation 
of state for the growing strings is incompatible with a constant dilaton, so the 
equations governing string cosmology. Equations (|^), effectively describe the 
spacetime. As approaches the string scale, the energy becomes extremely 
large; the expanding directions feel a large negative pressure and the contracting 
directions feel a large positive pressure, and the spacetime relaxes to a non- 
singular evolution. 

What in fact happens will depend on the correct expressions for p and p, 
the initial value of 0, and higher order corrections to the /3-function equations. 
Based on some simple computer simulations and qualitative arguments, we have 
found several possible outcomes for evolution via Equations (^. There seems to 
be little possibility for the rapid inflation and contraction to turn around before 
the expansion rate reaches unity in string units. What happens when it reaches 
this point depends on how we interpret the diverging integral in Equation (p6|). 
If we permit ourselves to imagine that the effect is a sufficiently large production 
of string matter then the expansion and contraction will suddenly reverse due 
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to the large pressures we expect to arise froni strings in rapidly expanding 
and contracting spacetimes. Once the turnaround occurs much will depend on 
the equation of state for less rapidly changing spacetimes. If expansion halts 
completely and the pressure drops to zero, then the spacetime will stabilize at a 
constant scale factor but the large energy density produced will cause the dilaton 
to grow to infinity in finite time, so that string loop effects become extremely 
important. We might imagine using the conjectured duality properties of string 
theory (see for example Sen 1994, Hull and Townsend 1995, Duff 1995, Witten 
1995) at this point to relate the background to the weakly coupled background of 
another theory. Another scenario occurs if the string coupling is too weak when 
the expansion rate becomes of order unity and not enough energy is produced; 
the large continued contraction of the angular directions can cause the otherwise 
growing dilaton to return to weak coupling as we reach the singularity. Of 
course, in the regime we are discussing all quantities - p, pk, \k, and so on - 
are of order one in string units; we can only speculate what the form of the 
/3-function equations and the equation of state might be. In discussing the form 
of p and p we should note also that at the point that the expansion becomes of 
order unity, the singularity in the Schwarzschild geometry is reached within a 
single string time. It is not clear that strings can equilibrate. Nonetheless, the 
blow-up of the integral in Equation (|9^) indicates that the description of the 
region near the putative singularity may be drastically modified by the quantum 
dynamics of string field theory. 

7 Conclusions 

We have argued that in regions of rapid spatial expansion, such as the region 
close to a spacelike singularity, a divergent energy density of string is produced 
as the curvature approaches the string scale. This leads us to one of two con- 
clusions, both of which are interesting and merit further investigation: either 
string theory the region close to the singularity undergoes a phase transition as 
the density of string approaches one string mass per string volume, or properly 
taking into account the back reaction of the produced strings on the spacetime 
pushes the singularity to infinite proper time. The answer could also be some 
combination of the two, or it could be that string theory is just inconsistent, 
although we feel that the divergence calculated has the flavor of a divergence 
arising from an incorrect approximation or extrapolation {i.e. from neglecting 
either backreaction or a phase transition at large densities), and so there is 
no obvious reason to call this divergence a symptom of some disease of the 
underlying theory. 

If the backreaction becomes important in the manner we have suggested, 
it is an indication that we must take seriously unstable macroscopic strings as 
sources in the /3-function equations. This not to say that solving the /3- function 
equations gives us an inconsistent picture, but rather that high order loop effects 
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become important in spacetimes with unstable strings. It would be interesting 
to find some sort of self-consistent way to include macroscopic strings as sources 
directly in the /^-function equation, so that one could evolve the spacctimc ac- 
cordingly and see what happens to singular regions and to the singularity the- 
orems. One might also envision trying to construct this approximation directly 
by resumming certain classes of loop corrections to the /3-function equations. 

Understanding how string theory modifies the dynamical equations of space- 
time is extremely important for a resolution of the problem of black hole evap- 
oration. A transition to a uniformly expanding region (Frolov et.al. 1990); con- 
tinued nonsingular contraction of the two-spheres and expansion of the radial 
direction (Martinec 1995) (where "radial" refers to the spacelike direction in the 
interior of the black hole); or simultaneous collapse of the expanding directions 
are all possible depending on the density of state and the equations of motion 
of the background. The first two possibilities tell us that the information is 
carried away to inaccessible regions of spacetime; the third possibility implies 
a slowly-evaporating-remnant scenario, where information is brought near the 
horizon by the recoUapse of the radial scale factor. 

String production in a non-trivial background spacetime should be an ap- 
proximation to the full quantum dynamics of the vacuum condensate. It may 
be that we have shown that it is incorrect not only in principle but in practice 
to separate dynamical quantum string fields and classical string backgrounds. 
If this is the case then we are stuck until we develop a manifestly background 
independent understanding of string field theory, or its nonperturbative com- 
pletion. 

Acknowledgements 

We would like to thank Eanna Flanagan, Vivek Iyer, and Martin O'Loughlin 
for enlightening conversations and helpful suggestions. A.L. and E.M. are sup- 
ported in part by funds provided by the DOE under grant No. DOE-FG02- 
90ER-40560. 

A Matching 4d black holes to 2d black holes 

Although we have argued that classical (tree level) string theory is insufficient 
for understanding spacelike singularities, it would be nice to better understand 
the classical geometry of singularities in string theory, in order to have some 
confidence in our starting point and to sharpen some of the arguments in the 
Introduction regarding the existence of singularities in string theory. We would 
like to repeat and extend some arguments made by Giddings et.al. (1993, 1994) 
concerning the worldshcct rcnormalization group flow of the a-model, here with- 
out the benefit of any gauge field or axion hair to stabilize the shrinking two- 
spheres. 
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Close to the singularity the spacctime metric should look like 

ds^ = -dT"^ + a^{T)dp^ + r^{T)dn'^ . (98) 

In the classical geometry wc approach the singularity in finite time as the metric 
in the p direction blows up, and we must go beyond the lowest order solution 
to the /3-function equations. It is well-known that the part of the a-model 
should flow to a trivial c = theory. Cecotti and Vafa (1992) have computed 
the Zamolodchikov metric along the entire renormalization group flow of the 
supersymmetric CP^ = S'^ model and found the metric to be well defined, 
even as becomes negative (see also Wittcn 1993). As the radius gets small 
and passes through zero the theory masses up and becomes trivial, with the 
soliton masses increasing as m = 8exp(— r^/2). Once the radius is sufficiently 
negative, the T, p plane will have excess central charge; a natural fixed point 
of the renormalization group flow is the SL{2,TZ)/U{1) black hole. This would 
mean that the velocity of the transverse stretching direction dominates the /?- 
function for a. One might convince oneself of this by adding a potential term to 
the spacetime action simulating the effects of string instantons wrapping around 
the 5^ direction (as do Giddings et.al. (1994)). This potential should be of the 
form 

V{R, (j)) = e-2'^-2-Va' (99) 

(Giddings et.al. 1994) where (p is the dilaton. This is a function of the scale 
factor which is of order unity at zero radius: it cannot compete with the large 

A terms present near a velocity-dominated singularity. 

This is an attractive scenario, but it has not been proven. Exact results 
such as those of Cecotti and Vafa (1992) do not help us much. The contribution 
of the 5^ degrees of freedom to the dilaton /3-function should be given by the 
Zamolodchikov c-function of the S'^ theory (Fradkin and Tscytlin 1985, Callan 
et.al. 1985, Polyakov 1986, 1987), which is only simply related to the Zamolod- 
chikov metric at RG fixed points (Zamolodchikov 1986). Even if wc had this 
information, cross-couplings between a and r would make the calculation ex- 
ceedingly difficult. 

We can also try to work backwards from the singularity, starting with the 
WZW black hole at small small r, (where r is the timelike direction inside the 
horizon, and t is the spacelike coordinate), times the S"^ theory close to its triv- 
ial infrared fixed point; with these initial conditions, the velocity terms coming 
from the coupling of the S'^ are small and we have some intuition for the evo- 
lution of the /3-function equations. Here we would use the following procedure. 
Begin with a spacelike slice S in this region. Given the mass of the WZW 
black hole, a starting point inside the horizon, and a starting point in the RG 
trajectory of the 5^ theory, we have a set of initial conditions for evolving the 
/3-function equations of the u-model away from the singularity in the radial di- 
rection. We may do this, for example, by using the well-known analogy between 
the /3-function equations in spacetime and the RG flow of the spatial part of the 
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a-model, where in the latter picture we would gravitationally dress the spatial 
c7-model with appropriate functions of the Liouville mode (Polchinski 1989, Das 
et.al. 1989, Banks and Lykken 1990). In the beginning of the radial cvolntion 
we find that the spacetime is some small perturbation of the WZW black hole 
coupled to the theory moving away from the trivial fixed point. If we begin 
with the S"^ theory close enough to this fixed point, as we pass through the 
horizon at r = r?^ the geometry will still be described by the WZW black hole 
times the massive theory with the mass slowly decreasing, and the r — t 
theory approaches the asymptotically fiat linear dilaton vacuum. At some point 
r = Tc the S'^ theory will begin crossing over to behavior characterized by the 
ultraviolet, asymptotically free fixed point and we should enter an asymptot- 
ically flat regime with a constant dilaton (since there is no longer any excess 
central charge to feed a linear dilaton solution). At this point, by Birkhoff's 
theorem, any deviations from the Minkowski metric should be that arising from 
the Schwarzschild solution at large radius. If we evolve the geometry from E 
towards the singularity, as long as the S"^ theory sufliciently close to its IR fixed 
point, we expect to reach the singularity of the WZW black hole, in finite time, 
since the effect of the flow of the S"^ theory should be minimal. The geometry 
is then that of a long throat with an effectively 2-dimensional geometry, termi- 
nating at one end at a 2d black hole singularity and opening out at the radius 
r = Tc, well outside the horizon of the 2d theory, into an asymptotically flat 
4-dimensional spacetime. It may be that when we reach Tc the crossover forces 
us into some other fixed point, but we believe this to be unlikely. 

Now let us move in the space of c = 4 conformal field theories by moving 
the initial position of the S"^ theory farther away from its trivial fixed point. 
This should have the effect of increasing the initial velocity of the coupling 
of the S"^ theory to r, which making cross-coupling terms between the radial 
and angular directions more important. We then expect Tc to approach rn] 
eventually along this marginal line, Th enters the asymptotically flat regime, 
unless we arc driven by the crossover behavior to some other fixed point. Since 
at this point, the velocity terms coming from at S shoiild be fairly large, 
wc no longer expect the region near the singularity to look like the WZW black 
hole: the RG flow might push the WZW singularity off to infinite distance, but 
we expect instead that we have some standard velocity- dominated singularity, 
based on our previous arguments. 
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